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ABSTRACT 

Context. Magnetic fields affect the local structure of the photosphere of stars. They can considerably influence the radiative properties 
near the optical surface, flow velocities, and the temperature and pressure profiles. This has an impact on observables such as limb 
darkening and spectral line profiles. 

Aims. We aim at understanding qualitatively the influence of small magnetic flux concentrations in unipolar plage regions on the 
centre-to-limb variation of the intensity and its contrast and on the shape of spectral line profiles in cool main-sequence stars. 
Methods. We analyse the bolometric and continuum intensity and its angular dependence of 24 radiative magnetohydrodynamic 
simulations of the near-surface layers of main-sequence stars with six different sets of stellar parameters (spectral types F to early M) 
and four different average magnetic field strengths (including the non-magnetic case). We also calculated disc-integrated profiles of 
three spectral lines. 

Results. The small magnetic flux concentrations formed in the magnetic runs of simulations have a considerable impact on the 
intensity and its centre-to-limb variation. In some cases, the difference in limb darkening between magnetic and non-magnetic runs is 
larger than the difference between the spectral types. Spectral lines are not only broadened owing to the Zeeman effect, but are also 
strongly affected by the modified thermodynamical structure and flow patterns. This indirect magnetic impact on the line profiles is 
often bigger than that of the Zeeman effect. 

Conclusions. The effects of the magnetic field on the radiation leaving the star can be considerable and is not restricted to spectral 
line broadening and polarisation by the Zeeman effect. The inhomogeneous structure of the magnetic field on small length scales 
and its impact on (and spatial correlation with) the local thermodynamical structure and the flow field near the surface influence the 
measurement of the global field properties and stellar parameters. These effects need to be taken into account in the interpretation of 
observations. 

Key words, magnetic fields - magnetohydrodynamics (MHD) - radiative transfer - stars: magnetic field - stars: low-mass - line: 
profiles 


1. Introduction 

Many cool mai n-sequence st ars have substantial surface mag¬ 
netic fields (see lRein^l20 1 2i for a review). These fields are be¬ 
lieved to be generated by a hydrodynamic dynamo in the deeper 
parts of the convection zone and are transport ed to the sur¬ 
face l ayers by buoyancy and convective flows (see ICharbonneaul 
1201 Ol for a review). Knowledge of the global geometry and 
strength of stellar surface magnetic fields and their dependence 
on stellar parameters is important in order to understand the un¬ 
derlying nature of th e dynamo pro cesses in stars of different in¬ 
ternal structure (e. g. lDonatil20i ll) . Moreover, the magnetic held 
is an important noise source for many astronomical obser vations 
(such as exo-planet detection, see e. g. lJeffers et alJ20l'3h . which 
makes a deeper understanding of its properties and the impact on 
observables essential. 

The presence of stellar magnetic fields is often indirectly de¬ 
tected by means of activity indicators (such as chromospheric S- 
index. X-ray flux) or photometric variability, which is mainly at¬ 


tributed to the rotational variation of the projected starspot area. 
Fast rotators in gen eral tend to be more magnetically active (see 
iReiners et al.ll201^ and references therein). Moreover, at least 
some stars were found to ex hibit activity cycles a nalogous to the 
11-year solar activity cycle (iBaliunas et al.llToO^ . 

In some cases, it is possible to measure the surface magnetic 
field by spectroscopic or spectropolarimetric means, which ex¬ 
ploit the action of the Zeeman effect on spectral line profiles 
in the integrated (Stokes-/) and circularly polarised (Stokes-V) 
light. While line broadening in Stokes-/ gives an estimate of 
the unsigned average magnetic held of a star, Stokes-E mea¬ 
surements are sensitive to the polarity of the held so that they 
pick up only the uncancelled net flux (signed average field). For 
main-sequence stars of spectral types F, G, and K, there exist 
only a few (Stokes-/) magnetic held measurements, as the mag¬ 
netic line broadening owing to the Zeeman effect is usual ly much 
small er than the rotational broadening in these stars (see lRein^ 
1201 21 and references therein). For M dwarfs, such measurements 
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are somewhat more feasible as the rotational velocities required 
for the generation of measurable magnetic fields are smaller. 
Some M dwarfs are found to be completely covered with a mag- 
netic field in the kG regime, with local values of up to ~ 7 kG 
(IShulvak et al.ll2014l) . 

For rapidly rotating stars, it is possible to infer the large- 
scale geometry of the magnetic field with spectropolarimetry: 
Zeeman-Doppler imaging (ZDI) exploits the latitude-dependent 
time dependence of the Zeeman-effect signature in polarised 
light caused by the Doppler effect due to the stellar rotation. ZDI 
reveals large-scale magnetic geometri es which are often quite 
unlike the solar one (iMorin et al.ll2010h . 

Owing to their interaction with the convective flows, stel¬ 
lar surface magnetic fields are expected to be highly structured 
and inhomogeneous on spatial scales as small as the scale of 
the convection, as it is known for the Sun: the magnetic flux 
at the solar surface is mainly concentrated in sunspots, pores, 
and smaller magnetic structures in the intergranular lanes. In 
these magnetic flux concentrations, the magnetic forces cause 
local modifications of the thermodynamical structure and the 
convective flows in the photospheric layers, which result, e. g. 
in a modulation of the local intensity emerging from magnetic 
flux concentrations. These modifications should be taken into 
account in spectroscopic and spectropolarimetric measurements 
of stellar magnetic fields as they will have a conside rable influ¬ 
ence o n the shape of the lo c al spe ctral line profiles. fBasri et aP 
(1 19901) and I Saar & Solankil (1 19921) analysed the significant im¬ 
pact of multicomponent atmospheres and magnetic field gradi¬ 
ents on spectral lines and the cons equences for stellar magnetic 
field measurements. More recently. iRosen & Kochukh^(l2012h 
showed that Zeeman-Doppler imaging of the magnetic field for 
a star with cool (dark) magnetic spots can fail if one does not 
simultaneously reconstruct the temperature. Three-dimensional 
MHD simulations of the near-surface convection as presented in 
this paper series provide a realistic and self-consistent structure 
of the magnetic field and stellar atmosphere and consequently 
can help to considerably improve stellar magnetic field measure- 
mentS; _ 

In iBeeck et alJ (1201 3ai hereafter Paper I), we presented non¬ 
magnetic 3D simulatio ns of six main-sequ ence stars (spectral 
types F3V to M2V). In iBeeck et al.l (12013bi hereafter Paper II), 
we analysed the impact of the convective flows and the sur¬ 
face structure caused by them on the limb darke ning and on 
the sh apes of the profiles of three spectral lines. In iBeeck et al.l 
(1201-51 hereafter Paper III) we investigated a grid of MHD sim¬ 
ulations with the same stellar parameters as the non-magnetic 
simulations analysed in Paper I and Paper II. The magnetic sim¬ 
ulation runs include a unipolar magnetic field with an average 
field strength of up to 500 G. In this paper, we analyse the im¬ 
pact of the magnetic field on the intensity field and spectral 
line profiles in the grid of simulations introduced in Paper I and 
Paper III. 

2. Methods 

The grid of six non-magnetic and 18 magnetic si mulations was 
run with the 3D r adiative RMHD code MURaM (IV6gleiil2003t 
IVogler et al.l 1200.51) . It comprises six sets of stellar parameters 
(Teff, log g), which roughly follow the (zero-age) main-sequence 
and correspond to stars of spectral types F3V, G2V, KOV, K5V, 
MOV, and M2V (see Paper I for more de tails). The metallic- 
ity wa s assumed to be solar (abundances bv I Anders & Grevessel 
Il989l) in all cases. The code use s the OPAL equation of state 
( Rogers et al.ll996HRogersl 19941) and ATLAS9 opacities (iKuruczl 


Il993l) . which were binned into four opacity bins (r-sorting ap¬ 
plying an averaged 3D atmosphere of the hydrodynamical run of 
each simulation as reference; for more details see Paper I). For 
each set of stellar parameters, there is one non-magnetic simu¬ 
lation and three magnetic simulations with average vertical field 
strengths of 20, 100, and 500 G (for more details, see Paper III). 
Table [T]lists the gravitational acceleration and the effective tem¬ 
perature of all 24 simulations considered. 

The spectral line synth esis was done with the forward mod - 
ule of the code Spinor (iFrutiged l2000t iFru tiger et al.l l200'oh . 
Whenever averaged spectral line profiles are considered, the av¬ 
erage is a spatial average over the 512 x 512 horizontal pixels 
(of the simulation boxes) as well as a temporal average over six 
snapshots which are taken several minutes apart from each other. 
These averaged line profiles were calculated for 10 different val¬ 
ues of /i (cosine of the angle between line of sight and stel¬ 
lar surface). For the disc-integrated line profiles considered, we 
used a numerical integration which takes into account differen¬ 
tial rotation according to 

fi(0) = Qeq (1 - a sin^ (1) 

where G is the local angular velocity of the star, Geq its value at 
the stellar equator, ff the latitude and a a parameter describing 
the differential rotation. 

In all cases considered here, the star was assumed to be ho¬ 
mogeneous on large scales as we use (six snapshots of) a local- 
box simulation for the line synthesis of the disc-integrated line 
profile. For more details, see Paper II and references therein. The 
simulation snapshots were taken with a cadence of a few minutes 
from the relaxed phase of the run (see Paper III). 

3. Intensity field 

3.1. Disc-centre intensity 

Figure [T] shows maps of the vertically emerging bolometric in¬ 
tensity from snapshots of the 500 G runs of the solar (G2V) and 
the M2V-star simulation. As shown in Paper III, the dark and 
bright structures in the intergranular lanes coincide with loca¬ 
tions of strong concentrations of magnetic flux, while the gran¬ 
ular upflows are only weakly magnetised. The high magnetic 
field strength in the flux concentrations causes a local depres¬ 
sion of the optical surface, which either results in bright struc¬ 
tures if the depth of the depression is comparable or larger than 
its horizontal extent (efficient side-wall heating) or dark struc¬ 
tures if its depth is small compared to the horizontal size. The 
smaller superadiabaticity in M dwarfs compared to hotter main- 
sequence stars is responsible for smaller typical depth-to-radius 
ratio of these local depressions. Consequently in F-, G-, and K- 
type dwarfs, there are only a few dark magnetic structures in 
the 500 G runs and none at all in the lOOG runs whereas the M 
star photospheres are dominated by dark magnetic structures at 
500 G and show many of them already at 100 G (for more details 
see Paper III). 

Table[T]lists the mean bolometric intensities at the disc cen¬ 
tre (i. e. yu = cos 7? = 1) along with its relative rms contrast and 
the extrema of a time-resolved sequence of intensity snapshots 
spanning several minutes for each of the 24 simulations consid¬ 
ered. The disc-centre intensity is almost unaffected by the mag¬ 
netic field in the F-, G-, and K-star simulations although the ra¬ 
diative flux (oc T*g) increases significantly owing to magnetic 
limb brightening (see Sect. 13.21) . In the M-star simulations, the 
disc-centre intensity is reduced in the 500 G runs with respect to 
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Table 1. Disc-centre (jj. - 1) intensities. 


Run (SpT - Bq) “ 

log.g[cgs] 

Teff [K] 

{^bol /^bol,0, hydro.) 

V<c>/<4oi> 

min(/boi)/(/boi)’’ 

max(/boi)/{/boi)’’ 

F3V - hydro." 

4.301 

6893 ± 6 

1.992 

0.206 

0.279 

2.031 

F3V-20G 

4.301 

6885 ± 6 

1.989 

0.203 

0.356 

2.238 

F3V- 100 G 

4.301 

6911+8 

1.994 

0.201 

0.392 

2.682 

F3V - 500 G 

4.301 

7003 ± 5 

2.017 

0.196 

0.345 

2.429 

G2V - hydro. >= 

4.438 

5764 ± 7 

1.000 

0.154 

0.533 

1.788 

G2V-20G 

4.438 

5779 + 9 

1.011 

0.152 

0.555 

2.368 

G2V- 100 G 

4.438 

5802 + 8 

1.015 

0.153 

0.555 

2.471 

G2V- 500G 

4.438 

5864 + 9 

1.009 

0.150 

0.425 

2.269 

KOV - hydro. >= 

4.609 

4856 + 6 

0.5055 

0.0796 

0.728 

1.840 

K0V-20G 

4.609 

4858 + 2 

0.5060 

0.0853 

0.739 

2.708 

KOV - 100 G 

4.609 

4878 + 4 

0.5125 

0.0945 

0.696 

2.632 

KOV- 500G 

4.609 

4901+2 

0.5099 

0.116 

0.505 

2.180 

K5V - hydro. 

4.699 

4368 + 2 

0.3235 

0.0705 

0.776 

1.687 

K5V-20G 

4.699 

4376 + 2 

0.3257 

0.0710 

0.802 

2.141 

K5V- 100 G 

4.699 

4383 + 3 

0.3270 

0.0734 

0.759 

2.128 

K5V - 500 G 

4.699 

4402 + 2 

0.3269 

0.0855 

0.624 

2.383 

MOV - hydro." 

4.826 

3905 + 1 

0.1949 

0.0359 

0.861 

1.203 

MOV-20G 

4.826 

3907 + 1 

0.1953 

0.0349 

0.865 

1.592 

MOV - 100 G 

4.826 

3909 + 1 

0.1955 

0.0350 

0.843 

1.580 

MOV- 500G 

4.826 

3906 + 1 

0.1929 

0.0538 

0.733 

1.433 

M2V - hydro." 

4.826 

3689 + 1 

0.1525 

0.0219 

0.904 

1.055 

M2V-20G 

4.826 

3691 + 1 

0.1526 

0.0218 

0.902 

1.283 

M2V- 100 G 

4.826 

3692 + 1 

0.1528 

0.0232 

0.865 

1.262 

M2V- 500G 

4.826 

3679 + 1 

0.1498 

0.0519 

0.787 

1.189 


“ SpT: spectral type; Bq: average vertical magnetic field strength 
global extrema occurring during a time-resolved sample of intensity snapshots spanning several minutes, 
non-magnetic (i. e. hydrodynamic) simulation. 
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Fig. 1. Vertically emerging bolometric intensity of a snapshot of the 500 G runs of the G2V {left) and M2V {right) simulations. The 
grey scale saturates at +2.5 standard deviations from the mean intensity. Intensity maps of further simulations are given in Paper I 
(Fig. 2) and Paper III (Fig. 1). 


the non-magnetic runs, whereas the radiative flux is unaffected 
(MOV) or only marginally reduced (M2V). 

In the F3V and G2V stars, the rms contrast is slightly re¬ 
duced hy the magnetic field while minimum and maximum in¬ 
tensity are somewhat increased as (moderately) bright magnetic 
regions replace the darkest regions in the intergranular lanes. In 
the K stars, the contrast is increasing with increasing Bq. Here, 
the bright magnetic regions are much more prominent than in 
the hotter stars. This is also reflected in the maximum value of 


the intensity, which is strongly enhanced in the magnetic K-star 
simulations to values up to 20 standard deviation above the mean 
intensity. The reason for this particular brightness of the bright 
magnetic regions is the weak dependence of the opacity on tem¬ 
perature between 4000 and 5000 K, which renders the location of 
the photospheric transition more sensitive to density than to tem¬ 
perature. Consequently, the partial evacuation in the magnetic 
flux concentrations shifts the optical surface to particularly high 
temperatures in these stars. A similar effect has been observed 
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Fig. 2. Intensity histograms of the G2V {left panels) and M2V {right panels) simulations. In both representations {upper panels: 
linear; lower panels: logarithmic), the histogram bins are linearly equidistant with a bin size of 0.01. 


in vertical vortices in the non-magnetic K-star simulations (see 
Paper II). From K5V to M2V, the bright magnetic regions be¬ 
come less bright (and less frequent) with decreasing effective 
temperature along our model sequence and the maximum inten¬ 
sity decreases (in terms of standard deviations above the mean). 

In the K-star and M-star simulations, the global intensity 
minima decrease from Bq = 20 G to Bq = 500 G, owing to the 
formation of dark magnetic regions. These dark regions become 
more pronounced with lower effective temperature as the depth 
of the optical surface depression becomes smaller compared to 
the horizontal size of the flux concentrations (see Paper III). 
Dark magnetic structures become dominating in the M-star sim¬ 
ulations and significantly reduce the mean disc-centre intensity 
in the 500 G runs with respect to the non-magnetic runs of these 
stars. 

Figure 12 shows histograms of the vertically emerging inten¬ 
sity (normalised to the mean) for all four runs of the G2 V and the 
M2V stars. The bright magnetic features produce high-intensity 
tails in the histograms in all magnetic runs of both stars. In the 
solar case, these bright features have the highest area fraction 
and intensity in the lOOG run and are considerably weaker at 
Bo = 20 G and Bq - 500 G. The dark magnetic structures are 
responsible for low-intensity tails in the histograms of the M2V- 
star simulations with Bq - lOOG and Bq = 500 G as well in the 
500 G run of the G2V star. As pointed out in Paper III, the ap¬ 
pearance of bright and dark magnetic structures does not only 
depend on spectral type, but also on the amount of magnetic 
flux available (and hence on Bq). Compared to the overall spread 
of the intensity distribution the dark structures are much more 


prominent in the M2V star simulations than in the G2V star (cf. 
Table[B. 

Figure[2shows a binned scatter plot of the vertical bolomet- 
ric intensity (normalised to its mean value) as a function of the 
magnetic held at the optical surface for all magnetic runs of the 
F3V, G2V, K5V, and M2V stars. In all runs (with the notable 
exception of the 500 G run of the M2V star), the mean intensity 
of the first bin (B(tr = 1) < 200 G), is close to the mean inten¬ 
sity of the whole box, while the intensity of regions with a sur¬ 
face magnetic held strength between 200 and lOOOG on average 
falls somewhat below this value. This has purely kinematic rea¬ 
sons: the magnetic flux is concentrated in the dark intergranular 
lanes, which stay dark as long as the held is too weak to con¬ 
siderably influence the thermodynamical structure and the flow 
patterns. At a higher held strength (B(tr = 1) > 1 kG), many of 
the curves have a positive slope, indicating that regions of high 
magnetic held strength have, on average, an increasing inten¬ 
sity with increasing field strength. This is the case for all runs 
for the F3V and G2V stars. In these two models, the 20 G runs 
show the highest intensities at a given held strength. The inten¬ 
sity of highly magnetised regions in the 500 G runs is somewhat 
lower than in the 20 G and 100 G runs and the scatter of intensity 
values (indicated by the error bars, which give the 1-cr scatter) 
within each bin increases somewhat. This can be attributed to 
the formation of a few dark magnetic structures among the still 
dominant component of bright structures. In the K5V simula¬ 
tion with Bq - 500G as well as in the lOOG and 500G runs of 
the M2V simulations, the slope of the binned scatter plots stays 
negative in the kG regime indicating that dark regions are dom¬ 
inant in these simulation runs. Only at very high held strengths 
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F3V G2V 





Fig. 3. Relation between the magnetic field strength B(tr = 1) at the optical surface and the vertical bolometric intensity /boi(/t = 1) 
normalised with its mean value. The field is binned in steps of 200 G. Error bars give the l-cr scatter within the bins. For bins with 
fewer than 1000 pixels (~ 0.06% of the total area considered), the symbol size is scaled with the logarithm of the number of pixels. 


of B(tr = 0) > 3 kG does the average intensity tend to increase 
with increasing field strength in these models. This is, however, 
produced by very few surface elements, which correspond to lo¬ 
calised bright points with very small area coverage. 

3.2. Centre-to-limb variation of intensity and contrast 

Figure |4] shows maps of the continuum intensity near 400 nm of 
the 500 G runs at an inclined viewing angle of jj - cos i? = 0.5. 
The snapshots are the same as in Fig. 1 of Paper III (and for G2V 
and M2V the same as in Fig. [T]in this paper). The differences 
between the different spectral types are mainly caused by the 
different depth of the depressions in the optical surface caused 
by the magnetic flux concentrations, but also by the different 
impact of the flux concentrations on the local density and tem¬ 
perature structures. In the F3V star, the depressions are deepest 
(in terms of pressure scale heights as well as in terms of gran¬ 
ule sizes, see Paper III) and the evacuation in the upper part of 
the magnetic flux concentration is strongest. Consequently, the 
intergranular lanes and the centres of the magnetic flux concen¬ 
trations are hidden from view behind granules while the granules 
which are seen through the highly evacuated upper parts of mag¬ 
netic structures appear brightened. Moreover, the granules show 
dark edges where the optical path through the cool layers above 
the optical surface is particularly long (“limb darkening” of sin¬ 
gle granules). With decreasing effective temperature the corru¬ 
gation of the optical surface as well as the typical depth of the 
depressions caused by the magnetic flux concentrations decrease 
(see Paper III) and the appearance of /(// = 0.5) changes. In the 


G-, K-, and M-star simulations, only a part of the granules seen 
through flux concentrations are brightened significantly. These 
sma ll scale brightening s are observed on the Sun as “faculae” 
(see iKeller et al.l12004 and references therein). In the M stars, 
especially M2V, the dark magnetic structures are still visible at 
ju = 0.5 and the faculae are less bright (in terms of the standard 
deviation of the intensity) than in the K- and G-type stars. 

Figure |5] shows the centre-to-limb variation of the contin¬ 
uum intensity in four narrow wavelength passbands for the G2V- 
and the M2V-star simulations, comparing the non-magnetic runs 
(dotted curves) with the 500 G runs. In both cases, the limb 
darkening is reduced in the magnetic runs compared to the 
non-magnetic runs, i. e. the magnetic field effectively brightens 
the limb. This limb brightening is produced by the faculae (cf. 
Fig. |4]i. The effect is stronger at higher effective temperature as 
the faculae become more prominent and have a higher filling 
factor. In the G2V-star simulation, the difference in limb darken¬ 
ing between the non-magnetic and the 500 G runs is even larger 
than the difference between the non-magnetic runs of different 
spectral types. 

In Figure |6] the centre-to-limb variation of the relative in¬ 
tensity fluctuations (“intensity contrast”) is plotted as a function 
of ju = cos 1 ? for the non-magnetic and 500 G runs of the so¬ 
lar (G2V) and M2V-star simulations. In the solar case, the in¬ 
tensity contrast at the disc centre does not significantly differ 
between non-magnetic case and 500 G case (cf. Table [T]i. Near 
the limb the effect of the faculae leads to a strongly enhanced 
intensity contrast: while in the non-magnetic case, the contrast 
monotonically drops from the disc centre to the limb, it reaches 
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Fig. 4. Continuum intensity in a passband between 400 and 410 nm of the six simulations with Bq - 500 G viewed at an angle of 
60° with respect to the vertical (i. e. ji - 0.5). 


a maximum in the 500 G runs around jj. - 0.3 in all wavelength 
passbands considered. In the M2V-star simulations, the magnetic 
field affects the intensity contrast at all disc positions: the con¬ 
trast is roughly twice (near the limb: three times) as high in the 
magnetic run as in the non-magnetic run. Unlike the solar case, 
the intensity contrast does not strongly increase from the disc 
centre towards the limb in the 500 G run as the faculae are less 
pronounced. 

Figure I?] shows the limb darkening in the non-magnetic and 
the 500 G runs of the G2V star in the Johnson UX, B, V, R, 
and I passbands (left panel) and the difference between the two 
runs in per cent (right panel). The highest relative deviations 
are in the UX band near the limb, where the magnetic run is 
more than 30% brighter than the non-magnetic run. For some 
real stars, limb darkening can be inferred from eclipses (in bi¬ 
naries) or exoplanet transits. Possibly, the limb darkening can 
be used as an indicator for small scale magnetic structures (fac¬ 
ulae) on such stars even in cases where the field is too small 
to be detected spectroscopically. In the same way as for the 
continuum passbands (see Fig. |3, the magnetic limb brighten¬ 
ing becomes weaker with decreasing effective temperature along 
our model sequence in the Johnson passbands (see Fig. lA.ll in 
AppendixlAli. 


4. Spectral line profiles 

4.1. Spectral line profiles at the disc centre 

In Paper II, we analysed the effects of the convective velocity 
field and its spatial correlation with the thermodynamical quan¬ 
tities on a set of three spectral lines (Fei at 616.5 nm, Fei at 
617.3nm, and Tit at 2223.3nm). Here, we investigate the im¬ 
pact of the local magnetic field structure and its correlation with 
thermodynamical quantities and with the convective flows on 
these lines. In the presence of a magnetic field, all three lines are 
split into their Zeeman components. The effective Lande factors 
are quite different: while the Fei line at 617.3nm was chosen 
for its large effective Lande factor of geff = 2.5, the Fe i line at 
616.5 nm is rather insensitive to the magnetic field (geit = 0.69). 
The third line, Tii at 2223.3 nm, has an effective Lande factor 
of .?eff = 1.66. However, as it is in the infrared, the broadening 
or splitting due to the Zeeman effect is expected to be larger in 
comparison to the line width than in the optical Fe i lines. 

Figure [8] shows a snapshot of the KOV-star simulation with 
Bo - 500 G. Four pixels are marked for which vertical local 
spectral line profiles of the Ti i line and the Fe i line at 617.3 nm 
are plotted in the two right panels. Point 1 is located in an upflow 
region with negligible magnetic field strength. The correspond- 
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G2V M2V 



Fig. 5. Limb darkening for the G2V (left) and M2V {right) simulations without magnetic held {diamonds and dotted curves) and with 
Bo = 500 G {squares and solid curves) in narrow passbands near 400, 600, 800, and 2000 nm. The temporal scatter is considerably 
smaller than the symbol size. 


G2V 


M2V 




Fig. 6. Centre-to-limb variation of the rms intensity contrast for the G2V {left) and M2V {right) simulations without magnetic held 
{diamonds and dotted curves) and with Bq - 500 G {squares and solid curves) in narrow passbands near 400,600, 800, and 2000 nm. 
The temporal scatter is considerably smaller than the symbol size. We note the logarithmic scale of the ordinates. 
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Fig. 7. Limb darkening in the Johnson UX, B, V, R, and I passband hlters (iJohnson & MorganI 195 It iBessellll 19901) in the G2V star. 
Analogous plots for the F3V and the MOV star are given in Fig. lA.ll in AnpendixlAl Left panel: Limb darkening of the non-magnetic 
{hydro.) and 500 G runs. Right panel: Relative difference between both runs. 


ing line prohles are shifted somewhat to the blue. Especially for 
the Fei line at 617.3 nm, this shift affects the line wings more 
strongly than the core, resulting in a strong asymmetry. As the 
wings of these photospheric lines form very close to the sur¬ 
face, while the line core forms somewhat higher in the stellar at¬ 


mosphere, this asymmetry reflects the deceleration of the (over¬ 
shooting) upflow as discussed in Paper II. Point 2 is located in 
a downflow with low held strength. The corresponding line pro¬ 
hles are shifted to the red, with an asymmetry (especially in the 
Fe I line at 617.3 nm) indicating an acceleration along the down- 
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2 3 4 

X [Mm] 


Fe I at 617.3 nm 



A-Xe [nm] 


Ti I at 2223 nm 



Fig. 8. Local vertical spectral line profiles for the KOV simulation with Bq - 500 G. Left panel: intensity image with four positions 
marked. Right panels: local profiles of the Fei line at 617.3 nm (top) and the Tit line at 2223 nm (bottom) at (p - 1); the four 
different colours correspond to the four points marked in the left panel. 


flow. We note that the continuum flux is lower in Point 2 and that 
the depth of the line profiles is smaller. In the average spectrum, 
Point 2 consequently has a lower weight than Point 1 for these 
lines. 

The other two points in Fig.[8]are situated in magnetised re¬ 
gions: Points 3 in a bright (hot) magnetic structure and Point 4 in 
the centre of a dark (cool) magnetic structure. The correspond¬ 
ing spectral line profiles are split into two components by virtue 
of the Zeeman effect. The components are asymmetric owing to 
the vertical gradient of the magnetic field: as the field rapidly de¬ 
creases with height in the photosphere (see Paper III), the cores 
of the Zeeman components are less strongly shifted away from 
the line centre than their wings. In comparison to Point 1, the 
equivalent width (EW) of the Fei line at 617.3 nm (upper right 
panel of Fig.|8ll is larger in Point 3 and smaller in Point 4, caused 
by the (moderate) temperature sensitivity of the corresponding 
line opacity. As the continuum is also higher in Point 3 than in 
Point 4, the weight of the local line profile from Point 3 will be 
higher in the spatially averaged line profile than that of Point 4. 

In contrast to the Fei line at 617.3nm, the Tii line (lower 
right panel in Fig. [8]l has a strongly reduced EW in the bright 
magnetic structure (Point 3) with respect to Points 1 and 2. The 
Ti I line is very sensitive to temperature around 4000 K, owing 
to ionisation and to depopulation of the lower level of the tran¬ 
sition. Therefore, the line opacity is considerably reduced in hot 
magnetic structures of this star. This effect will be referred to 
as “line weakening” in the following discussion. The same phe¬ 
nomenon also exists in the Ee i lines but at a higher temperature 
(around 6000 K), which is more relevant for the G2V- and E3V- 
star simulations. In the dark magnetic structure (Point 4), the 
EW of the Ti i line is somewhat larger than in Points 1 and 2. In 
summary, dark (cool) and bright (hot) magnetic structures can 
produce crucially different local spectral line profiles. As shown 
in Sect. 13.11 the relative importance of dark and bright magnetic 
structures depends on the average field strength and the stellar 
parameters. 


Pigure0 shows horizontally averaged profiles of the Ee i line 
at 617.3 nm for the 500 G runs of the E3 V- and the KOV-star sim¬ 
ulations and of the Ti i line for the 500 G runs of the KOV- and 
the M2V-star simulations. The thick black curves represent the 
line profiles averaged over the entire area of the simulation box, 
while the four coloured curves are the profiles averaged over 
surface components with different magnetic field ranges at the 
optical surface. The numbers in parentheses give the weight of 
the individual components in the average line profile. The pro¬ 
file of the Ee i line at 617.3 nm has a strongly extended red wing 
in the E3V simulations. As the decomposition into distinct sur¬ 
face components illustrates, this red wing is produced mainly 
by the flow velocities in the area components with intermedi¬ 
ate field strengths of 250 to 2000 G: these components mainly 
represent short-lived, dynamic magnetic flux concentrations in a 
convective collapse phase in this star (see Paper III). This results 
in an (on average) increased downflow velocity in the optical 
depth range where this spectral line is forming. Additionally, 
the Zeeman effect broadens the line profiles emerging from 
these surface components, shifting its red wing even further to¬ 
wards longer wavelengths. The regions with high magnetic field 
strengths of > 2 kG do not contribute much to the average spec¬ 
trum because of line weakening, which affects the Ee i lines in 
the G2V and E3V atmospheres (while it does not play a great 
role for KOV, cf.Eig.|8]l. 

In the KOV simulation, the flow velocities are much 
lower and the magnetic flux concentrations are less dynamic. 
Consequently, the Doppler shift of the line components is 
weaker. As the temperature is lower, less line weakening occurs 
in the Ee i line at 617.3 nm and the substantially broadened line 
profile emerging from the strongly magnetised area components 
produces extended wings in the averaged (disc-centre) line pro¬ 
file. Eor the Ti i line, these wings are even more pronounced as 
the Zeeman splitting is stronger in the infrared (with respect to 
the line width). Owing to line weakening, however, the equiv¬ 
alent width of the line profiles stemming from the magnetised 
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F3V - Bo = 500G 


Fe I (617.3 nm) 


KOV 



KOV - Bq = 500G - Ti I (2223 nm) 



Bo = 500G 


Fe I (617.3 nm) 



M2V - Bo = 500G - Ti I (2223 nm) 
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Fig. 9. Decomposition of four line profiles (at /r = 1) into differently magnetised components. The black curve shows the average 
line profile, while the coloured curves show the line averaged over areas where B(tr = 1) is in the specified ranges. Each line is 
normalised to its own continuum. The numbers in parentheses give the weight of the components in the composite spectrum (i. e. 
average continuum flux times area fraction). 


F3V - Fe I at 617.3 nm F3V - Fe I at 617.3 nm 




F3V - Fe I at 617.3 nm F3V - Fe I at 617.3 nm 




Fig. 10. Centre-to-limb variation of the equivalent width, FWHM, and the effective Doppler shifts of the line core and wings (from 
left to right) of the Fe i line at 617.3 nm in all four simulation runs of the F3V star. The colour code is the same for all panels, see 
legend in the left panel. 


area is somewhat reduced. In the M2V simulation, line weaken¬ 
ing does not play a role for the lines considered. The line width is 
larger here for all surface components owing to the higher pho- 
tospheric pressure, rendering the extended wings more difficult 
to see in the pressure-broadened line profile. 

4.2. Centre-to-limb-variation of spectral line profiles 

In Paper II, the variation of the equivalent width (EW), the full 
width at half maximum (FWHM), and the Doppler shifts of the 
cores and wings of spectral lines (in particular of the Fe i line at 
617.3 nm) as a function of viewing angle p - cost? were anal¬ 
ysed for the non-magnetic simulations. Here, we show the ef¬ 


fect of the magnetic field on the centre-to-limb variation of the 
line profiles. As discussed above for the vertical (i. e. disc-centre) 
spectral line profiles, the magnetic field has a substantial impact 
on the line profile shape parameters. Figure[T0]shows this for the 
Feiline at 617.3 nm for the example of the F3V star. The EW de¬ 
creases with increasing Bq at all disc positions (i. e. for all values 
of p). This is mostly caused by line weakening in the magnetic 
flux concentrations (cf. Fig. [8]). The centre-to-limb variation of 
the FWHM does not differ much between the non-magnetic and 
the 20 G and 100 G runs. At Bq — 500 G, however, the FWHM 
near the limb is considerably smaller than in the non-magnetic 
run. This is a result of the stronger corrugation of the optical 
surface in the magnetic run, which hides from view the part of 
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Fig. 11. Difference of equivalent width, FWHM, line core shift, and line wing shift (from left to right) between simulations with 
Bq - 500 G and the respective non-magnetic simulations for the three spectral lines as a function of /i. 


the granular flows receding from the observer as well as the cool 
downflows and thus reduces the broadening induced by varia¬ 
tions of the line-of-sight velocity and temperature. The centre- 
to-limb variations of the effective Doppler shifts of the line dif¬ 
fer considerably between the runs. At the disc centre, the line 
wing^ are strongly shifted to the red in the 500 G runs owing 
to the higher downflow velocities at constant optical depth (cf. 
Fig|9]l in this run in comparison to the non-magnetic runs. Near 
the limb, however, the stronger corrugation in the 500 G run re¬ 
sults in an effective blue shift, as the receding flows are mostly 
hidden behind the granules. As the effect of the magnetic field 
on the shifts of line cores and wings is qualitatively and quanti¬ 
tatively different, the magnetic field also has a strong impact on 
the line bisectors. 

* Consistent with the definition in Paper II, the line wing shift is de¬ 
termined as the displacement of the bisector from the rest wavelength 
of the spectral line at 5% of the line depth. 


Figure [TT] shows the difference in EW, FWHM, line core 
shift, and line wing shift between the 500 G runs and the non¬ 
magnetic runs as a function of p for all spectral types and all 
three spectral lines (for the Fei line at 617.3 nm in the F3V 
star this corresponds to the difference between the blue and the 
black curve in Fig. [TOl i. The differences in EW and EWHM are 
given as relative difference in per cent, while the Doppler shifts 
are given as absolute differences in km s '. Eor both Eei lines, 
the EW is smaller with magnetic field than without in the G2V 
and E3V simulations (cf. Eig. fTOl i owing to line weakening in 
the bright magnetic flux concentrations. In the other stars, line 
weakening does not play an important role for these lines. Here, 
the EW is unaffected in the magnetically insensitive Ee i line at 
616.5 nm while it is increased in the magnetically sensitive line 
at 617.3 nm owing to broadening by the Zeeman effect, which 
reduces the degree of saturation of the line. Eor the Ti i line, line 
weakening also affects the EW in the KOV star at disc centre. 
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F3V - Fe I at 616.5 nm F3V - Fe I at 617.3 nm F3V - Ti I at 2223.3 nm 
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Fig. 12. Doppler shift of the line profile cores (solid curves) and wings (dashed curves) of the three spectral lines in the F3V 
simulations as function of itrot sin i. In all cases, i - 60° and a - 0.2 was assumed. 


We note that the EW of this line is very small in the F3V sim¬ 
ulation, particularly near the disc centre (see Paper II) and only 
regions with exceptional temperature gradients and high density 
contribute substantially to the average profile, which might ex¬ 
plain why the impact of the magnetic field on the EW of this line 
is not similar or stronger than in the G2V star as expected. 

The EWHM follows a similar trend; for the hotter stars, the 
EWHM tends to be smaller with magnetic field while, in the 
cooler stars, it is almost unchanged in the magnetically insen¬ 
sitive line and strongly enhanced in the magnetically sensitive 
lines. This again is showing the relative importance of line weak¬ 
ening and Zeeman broadening. As the line weakening affects 
mainly the contributions from downflows, it reduces the EWHM. 

The effective Doppler shifts in the lines show an impact of 
the magnetic field which is almost independent of the sensitivity 
of the line to the magnetic field. The disc centres of the F3V and 
G2V stars show strongly redshifted line wings and slightly red- 
shifted line cores in the 500 G runs (relative to the non-magnetic 
runs) in both Fei lines. As shown earlier, this is caused by a 
strongly enhanced red wing, resulting from the larger filling 
factor of downflows (see Paper III) and the higher downflow 
speed at the optical depth range where these lines form. The 
Zeeman effect only slightly enhances this effect by broaden¬ 
ing the spectral line originating from the magnetised downflows. 
Consequently, the relative redshifts in line wings and cores are 
somewhat higher in the 617.3 nm line than in the 616.5 nm. 
Towards the limb, these two stars show a small relative blueshift 
of line wings and cores with respect to the non-magnetic runs 
(which is probably due to the stronger corrugation of the opti¬ 
cal surface). In the cooler stars, the impact of the magnetic field 
on the Doppler shifts of the line cores and wings and on their 
dependence on p is much weaker since the convective velocities 
(particularly in the magnetic flux concentrations) are consider¬ 
ably smaller. 

4.3. Disc-integrated spectral line profiles 

Applying the numerical method outlined in Paper II, disc- 
integrated spectral line profiles were calculated using the an¬ 


gle dependent spectral line profiles from the local simulation 
boxes. So far, we do not include any large-scale surface inho¬ 
mogeneities (such as active regions), but assume the stars to be 
homogeneous on large scales (i. e. having the properties of the 
3D simulation with a constant Bq across the entire disc). The 
disc integration includes the rotation of the star, which follows 
the simple law given in Eqn. ([T) and thus introduces three param¬ 
eters: the equatorial velocity of the rotation, Vrot, the inclination, 
/, of the rotation axis with respect to the line of sight, and the 
differential rotation parameter, a. 

As shown in Figs. [10] and (TT] the 500 G runs of the hottest 
two stars of our sequence show redshifted line wings and (to 
a lesser extent) redshifted line cores at the disc centre. In the 
disc-integrated spectral line profiles, there is still a relative red- 
shift with respect to the non-magnetic simulations, but the shift 
is considerably smaller. These effective Doppler shifts turn out 
to be sensitive to the rotational velocity as the spectral line pro¬ 
files get distorted by the rotational broadening (cf. Paper II). 
Figure [12] illustrates this for the F3V star (where the shifts are 
largest) with the assumption of i = 60° and a - 0.2. In all three 
spectral lines, the line core moves redwards with increasing rota¬ 
tional velocity and the line wings move bluewards. Although this 
is the case in the non-magnetic and the magnetic runs alike, the 
effect becomes stronger with magnetic field; in the 500 G runs, 
the redshift in the line core reaches values around 1.5 km s * at 
a L'rot sin i of 15 km s“^ 

From the above (Sects. I4~n and l4.2l i it is evident that the mag¬ 
netic field has an impact on the spectral line profiles not only by 
virtue of the Zeeman effect, but also owing to the local modifi¬ 
cations of the thermodynamical structure and the flow patterns 
caused by the magnetic flux concentrations. In order to disen¬ 
tangle these two effects, the spectral line synthesis for the 500 G 
runs was carried out again ignoring the magnetic field for the line 
synthesis. The result is indicated by the red curve in Fig.[T3]for 
the example of the Fe i line at 617.3 nm and the G2V-star simula¬ 
tion. For comparison, the line profile of the non-magnetic simu¬ 
lation (grey) and of the full synthesis of the 500 G run (blue) are 
shown as well. Without rotation, the difference between the grey 
and red curves, which is exclusively caused by the modifications 
in thermodynamical structure and flows, is approximately equal 
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G2V - Fe I (617.3 nm) - no rotation 



G2V - Fe I (617.3 nm) - sin i = 7.5 km s" 



Fig. 13. Disc-integrated profiles of the Fei line at 617.3nm without rotation {left) and with u^ot sin/ = 7.5kms ' (differential, 
a = 0.2, i - 60°; right) for G2V simulations. The grey line represents the profile of the non-magnetic run {hydro.), the blue curve is 
the profile of the 500 G run, the red curve is calculated using the thermodynamic structure of the 500 G run but assuming B = QG for 
the line calculation (see text for further explanation). Dashed curves give the bisectors of the profiles with A - values magnified 
by a factor of 5 for clarity. We note the different scales of ordinate and abscissa. 


in magnitude to the change induced by the Zeeman effect (dif¬ 
ference between blue and red curves). While the local structure 
modifications shift the line somewhat to the red and reduce its 
equivalent width and FWHM, as discussed above (cf. FigJTTli. 
the Zeeman effect broadens the line, as expected. This broad¬ 
ening, however, is not exactly symmetric. The magnetic field is 
predominantly located in downflows and is thus correlated with 
a Doppler shift leading to a more broadened red wing. As a re¬ 
sult the Zeeman effect increases the C-shape of the profile of this 
spectral line as can be seen in the bisectors (dashed curves in the 
figure). In the other stars and in the Ti i line a qualitatively sim¬ 
ilar impact on the bisectors is visible (see Figs. IA.2l and IA.3b . 
albeit less pronounced. 

At a moderate rotational velocity of 7.5kms“' (see right 
panel of Fig. [T3l l. the broadening due to the Zeeman effect 
is already very small compared to the rotational broadening. 
Consequently, the effect of the modified flows and thermody¬ 
namical structure on the line is overall much larger than the 
Zeeman effect. 

This calculation has been carried out for all simulations and 
all spectral lines (more plots are given in Appendix lAli. We find 
the locally modified thermodynamics and flows to have a major 
impact on the line profiles in the F-, G-, and K-type stars. Only 
in the M-star simulation is the effect of thermodynamics signif¬ 
icantly smaller than that of the Zeeman effect. This has several 
reasons. Firstly, line weakening does not play a role here for the 
spectral lines considered. Secondly, the absolute contrast of the 
temperature is smaller than in the hotter stars. Thirdly, the con¬ 
vective velocities are generally smaller; in particular, there are 
no large downflow velocities in the magnetic structures as they 
are found in the hotter stars. 

5. Concluding remarks 

We have shown that the intensity, its contrast and centre-to-limb 
variation as well as spectral line profile shapes are strongly in¬ 
fluenced be a moderate local magnetic field. In unipolar regions 
with an average magnetic field of up to 500 G, the magnetic flux 
is strongly concentrated in localised structures of reduced or en¬ 
hanced intensity. F-, G-, and K-type stars show more magnetic 
flux concentrations of enhanced intensity under these condi¬ 
tions, whereas M-type stars have prominent dark pore-like struc¬ 
tures. These structures generally become bright near the limb 
and thus reduce the limb darkening (magnetic limb brighten¬ 
ing). Moreover, the intensity contrast is enhanced, particularly 


near the limb. The influence on spectral line profiles is not only 
caused by the action of the Zeeman effect on the magnetically 
sensitive lines, but also by the presence of the magnetic flux con¬ 
centrations which substantially deviate in their thermodynamical 
structure and flow patterns from the non-magnetised surround¬ 
ings. Moreover, the height dependence of the magnetic field and 
the flow velocity as well as convective collapse events (resulting 
in strong downflows in magnetic flux concentrations) add asym¬ 
metries to the full-disc spectral line profiles, which can strongly 
depend on the rotation velocity of star. 

For the synthesis of spectral lines as well as for the calcu¬ 
lation of the limb darkening, we have implicitly made the as¬ 
sumption that the stars are completely covered with a unipolar 
field, which is homogeneous on large scalesQ As the solar ex¬ 
ample as well as ZDI results show, stellar magnetic fields are not 
only structured on small scales, but also possess global structure 
(e. g. active latitudes, star spot groups, polar fields). Even though 
a homogeneously magnetised star as we implicitly assume does 
probably not exist in reality, the results presented in this paper 
illustrate that the impact even of a moderate magnetic field on 
observable quantities is expected to be significant. For Stokes-/ 
measurements the Zeeman-effect signature can be largely oblit¬ 
erated by the effect of the deviating atmospheric structure in the 
magnetised regions. Consequently, magnetic field measurements 
for stars (at least for solar-like stars) have to be taken with cau¬ 
tion because the magnetic field is usually considered to have no 
impact on the spectral lines apart from the Zeeman effect. The 
com parison of magn etically sensitive and insensitive lines (see 
e. g. lRobins(Mll98flh only helps to a certain extent, as this does 
not account for the correlation of the magnetic field with, e. g. 
the velocity field, which can modify the impact of the Zeeman 
effect. 

In this paper, we do not consider spectral line profiles in po¬ 
larised light. The line weakening discussed here does also affect 
the other Stokes components of the spectral lines. Moreover, a 
strong line-of-sight velocity variation in the magnetised surface 
components (as we have found it at least for the F3V star) will 
also affect the Stokes-V signal and reduce the resolution of ZDI 
(which makes use of the line-of-sight velocity to obtain informa¬ 
tion on the spatial distribution of the magnetic flux). 

Great progress has been made in the detection and measure¬ 
ment of stellar magnetic fields in recent years. To improve the 


^ As we only consider the Stokes-/ component, it is not an important 
issue that such a global configuration is not divergence-free. 
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methods it will be important to include the various effects of the 
magnetic field on the stellar atmospheres resulting from its 3D 
structure. 
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impact on the spectral line profile. In the K stars, the thermody¬ 
namic modifications seem less important than the Zeeman effect 
without rotation (left panels). In the rotating case (right panel), 
both effects are of similar magnitude. As they are of opposite 
sign, the apparent broadening of the line by the Zeeman effect is 
reduced by the modifications in the atmosphere structure. In the 
M stars (here, only the MOV is shown), the impact of the mod¬ 
ified thermodynamics and flows is very small. For these stars, 
systematic errors in Stokes-/ measurements of the magnetic field 
caused by differences in atmospheric structure between magne¬ 
tised and unmagnetised regions will probably be the small. 

Figure |A3] shows the same plots as Fig. lA.2l but for the Tii 
line and only for the K- and M-type stars. Although the line 
weakening is stronger for this line in the K-type stars (cf. Fig. Ell, 
possibly leading to an underestimation of the unsigned average 
field from Stokes-/ measurements, the general detection of mag¬ 
netic field is probably more feasible in this line as the broaden¬ 
ing due to the Zeeman effect is larger and leads to extended line 
wings. The line profiles from the M-star simulations are almost 
unaffected by the modifications of the thermodynamical atmo¬ 
sphere structure. 


Appendix A: Additional figures 

In Fig. |7] in Sect. 13.21 the magnetic limb brightening is given 
in the Johnson passbands UX, B, V, R, and I for the G2V star. 
Figure lA.ll shows the limb brightening for the F3V and MOV 
simulations. The magnetic limb brightening is largest for the 
F3V star and becomes smaller towards lower effective temper¬ 
atures as the faculae become less bright and less frequent (cf. 
Fig.lUi. 

In Fig.[T3in Sect. 14.31 disc-integrated profiles of the Fe i line 
at 617.3 nm are given for the G2V simulations without magnetic 
field (grey curves) and with Bq = 500 G (blue and red curves) 
where the magnetic field has been ignored for the spectral line 
synthesis for the red curves in order to disentangle the effect 
of the modified atmospheric structure from that of the Zeeman 
effect. Here, we give analogous plots for further simulations for 
the same spectral line (Fig. IA.2I) and for the Ti i line at 2223 nm 
(Fig.lAj. 

As Fig. lA.21 illustrates, the Fei line is significantly affected 
by line weakening in the F3V star. Similar to the G2V-star case, 
the line is shifted to the red (in particular its wings) by the mod¬ 
ified convective flows. The Zeeman effect has a relatively small 


13 


I(m)/I(1) I(m)/I(1) 


B. Beeck et al.: 3D simulations of near-surface convection in main-sequence stars. IV. 


F3V 



MOV 



F3V 



MOV 



Fig.A.l. Same as Fig.|7]but for the F3V and the MOV star. 
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B. Beeck et al.: 3D simulations of near-surface convection in main-sequence stars. IV. 


F3V - Fe I (617.3 nm) - no rotation F3V - Fe I (617.3 nm) - sin i = 7.5 km s 




KOV - Fe I (617.3 nm) — no rotation KOV - Fe I (617.3 nm) - sin i = 7.5 km s 




K5V - Fe I (617.3 nm) — no rotation K5V - Fe I (617.3 nm) - sin i = 7.5 km s 




MOV - Fe I (617.3 nm) - no rotation MOV - Fe I (617.3 nm) - sin i = 7.5 km s 




Fig. A.2. Same as Fig.[T3l but for the F3V, KOV, K5V and MOV simulations. 
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B. Beeck et al.: 3D simulations of near-surface convection in main-sequence stars. IV. 


KOV - Ti I (2223 nm) - no rotation 



K5V - Ti I (2223 nm) - no rotation 



MOV - Ti I (2223 nm) - no rotation 



M2V - Ti I (2223 nm) - no rotation 



KOV - Ti I (2223 nm) - sin i = 7.5 km s" 






Fig. A.3. Same as Fig. lA.2l but for the Ti i line at 2223 nm and for the K- and M-star simulations. 
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